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Abstract 

Here we investigate the holographic dark energy 
model in the framework of FRW cosmology where the 
newtonian gravitational constant, G, is varying with 
cosmic time. Using the complementary astronomical 
data which support the time dependency of G, the evo- 
lutionary treatment of EoS parameter and energy den- 
sity of dark energy model are calculated in the pres- 
ence of time variation of G. It has been shown that 
in this case, the phantom regime can be achieved at 
the present time. We also calculate the evolution of G- 
corrected deceleration parameter for holographic dark 
energy model and show that the dependency of G on 
the comic time can influence on the transition epoch 
from decelerated expansion to the accelerated phase. 
Finally we perform the statefinder analysis for G- cor- 
rected holographic model and show that this model has 
a shorter distance from the observational point in s — r 
plane compare with original holographic dark energy 
model. 



keyword: Cosmology, dark energy, holographic 
model, gravitational constant. 



1 Introduction 

Nowadays, it is widely believed that the cosmos is ex- 
periencing an accelerated expansion. This idea and 
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belief came into existence after collection of data 



from "Type la supernova" in 1998 (iPerlmutter et al 



1998). Also the other data from WMAP (Bennett et al 



2008 ),SDSS dTegmark et al Il2004h and X-rav dAUen et al | 



20041 ) experiments support this accelerated expansion. 



In the framework of standard cosmology, the existence 
of dark energy with negative pressures is essential to 
interpret the cosmic acceleration. Hence, dark energy 
scenario has got a lot of attention in modern cosmology 
both from theoretical and observational point of view. 
Observationally, The result of SNela experiment shows 
that dark energy occupies about 72% of the total energy 
of our universe, dark matter and ba ryons about 28% 



of th e total energy of the universe (jPerlmutter et al 



19981) . Although the nature of dark energy is still 
un-known, but the ultimate fate of the current uni- 
verse is determined by this mysterious component. Till 
now, some theoretical models have been proposed to 
interpret the behavior of the dark energy.The first and 
the simplest one is Einstein 's " cosmological constant" 
(jSahni fc Starobinskv Il2003( ) which, of course, has two 
problems called fine -tuning and cosmic coincidence. 
The cosmological constant has the fixed equation of 
state WA = — 1, while the dynamics of current expan- 
sion can be explained by dynamical dark energy mod- 
els with time varying equat ion of state . The scalar 



fields such as quint essence (jWetterich I I1998I ). phan 



tom ICaldweU hOQ^ or the combination of both which 
is called quintom ( Elizalde et al 20041 ) are examples 
of dynamical models. The other dynamical dark en- 
ergy models which interprets the current acceleration 
of expansi on are constru cted based on quantum grav- 
ity theorv dWitten 1 120021 ). The holographic dark en- 



ergy (HDE) model is one of the the proposed models 
base d on the holographic p rinciple in quantum grav- 



er 

ity (jHorava fc Minic 1 12000) ■ According to the holo- 
graphic principle, a short distance ultra-violet (UV) 
cut-off is related to the long distance infra-red (IR) 
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cut-o ff, due to the limit set b y the formation of a black 



hole iHorava fc Minic I boOOl) . The holographic prin- 



ciple indicates that the zero-point energy of a system 
with size L should not be exceed from the mass of black 
hole with the same size. From the above principle, the 
energy density of HDE model in cosmology can be de- 
scribed as: 



3c^ 



Pd 



(1) 



Where L is the cosmic horizon and c is a numerical 
constant of order unity and G is a Newton's gravita- 
tional constant. The length scale L has an essential 
role in the definition of energy density of HDE model. 
Therefore the various model of HDE have been con- 
structed for different of infrared (IR) cutoff length. For 
example the simple choice of IR cutoff is the Hubble 
length which leads to w rong equation of state for DE 
( Horava fc Minic IboOO ) . However in the presence of in- 
teraction between dark matter and DE the HDE model 
with Hubble radios IR cutoff can derive the acceler- 
ated expansion and also solv e the coincidence prob- 



lem (jPavon &: Zimdahl 1 120051 ) . The other choice for 
IR cutoff is the particle horizon. In this case, like 
Hubble l ength, the accelerated ex pansion cannot be 



achieved iPavon fc Zimdahl I ( 20051) . Another choice is 



the event horizon where the cosmic acceleration can be 



interpreted in this case (jZhou et al 1 120071 ) . Nojiri and 
Odintsov (2006) investigated the holographic DE model 
by assuming IR cutoff depends on the Hubble rate, par- 
ticle and future horizons, span of life of the universe 
and cosmological. In this generalized form of HDE the 
phantom regime can be achieved and also the coin- 
cidence problem is demonstrated. Unification of early 
phantom inflation and late time acceleration of the uni- 
verse is the other feature of this model. Recently, the 
HDE model has been constrained by var i ous astronom 



ical observations (|Hua.ng fc Gong 



20051 : IWu et al. Il2008t Enqvist et al 



20041: IZhang fc Wu 
20051) . 



In addition, there are some theoretical and obser- 
vational supports indicating that Newton's gravity 
constant varies and changes with cosmic time. The 
first theoretical i dea in this r espect is the pioneer- 
ing wo rk of E)irac(pirac 19381) . and then the idea of 
Dvsonl Dvson I Il972l) . Also, the Branse-Dicke frame- 



work in the Physics predicts t he va riability and fiuc- 



tuation of G ljBrans fc Dicke I 119611 ). Moreover the 



varying behavior of G in Kaluza-Klein theory was 
associated with a scalar field appearing in the met- 
ric component corres ponding to the 5-t h dimension 
and its size variation (jKaluza. et al 1 119211 ). In this the- 
ory, a scalar field paired with gravity by a new pa- 
rameter replaces the quantum gravity. The variabil- 



ity of G with time is also supported from the obser- 
vational viewpoint. The observational data collected 



by Type la S upernova (IGa ztanag a et al 1 120021 ) .Hulse- 
Taylor Binarv (|Damour et al ■■1988,1 . astro-s eismological 



data from pulsating white dwarf stars ( Benvenuto 



2004 : lBiesiada fc Malec 12004 ). helio-sesmiolo gical ( Guenther"^ 
1998) and the Big Bang Nuclei-synthesis datal Copi et al ^ 
20041) support a variable value for G with time. We re- 
fer to these observations in the section 3 of the paper. 
Here in this work we consider the HDE model with 
time varying G, the so-called G-corrected HDE model, 
in spatially flat FRW universe. We consider the event 
horizon as an IR cut-off in relation ([1]). In this concern 
some other works have been done in which the HDE 
model h as been considere d with time dependency of 
G, i.e., ( Jiano et al 20091) . Here we obtain the equa- 



tion of state Wd as well as deceleration parameter q and 
stateflnder pair {s, r} for G-corrected HDE model in 
FRW universe and also solve the related equations nu- 
merically by using the observational values for G{t). 
It is clear that constraining a given model against the 
observational data is model dependent. Therefore some 
doubts usually remain on the validity of the constraints 
on cosmological parameters. In order to solve this prob- 
lem, we use the cosmography, i.e. the expansion of 
scale factor in Taylor scries with respect to the cos- 
mic time. The flrst term of Taylor series is the Hubble 
parameter (H = ^), the second term is the decel- 

,2 

eration parameter {q = —^^tts^^j), the third term is 
the jerk parameter [r = aH'^'dt-^ )' forth term is 

,4 

snap (fc = ^^4°^j ) and the flfth term is Icrk param- 
eter {I ~ aH'-dt'^ )• "^^^ present values of the above 
parameters can be used to describe the evolution of 
the universe. For example qo < indicates the cur- 
rent accelerated expansion of the universe and also ro 
allows to discriminate between different dark energy 
models. Since Bubble's parameter which corresponds 
to the first derivative of the scale factor (a) and the de- 
celeration parameter which corresponds to the second 
derivative of the scale factor (a) can not distinguish 
between the different models, we need a higher der iva- 
tive of scal e factor. Sahni et a l. (|Sahni et al Il2003l) and 
Alam et al.( Alam et al 2003b[ ). by using the third time 
derivative of scale factor, introduced the statefinder 
pair {s, r} in order to diagnosis the treatment of dark 
energy models. The statefinder pair in spatially fiat 
universe is given by: 



r-1 



(2) 



The statefinder parameters s and r are the geometrical 
parameters, because they only depend on the scale fac- 
tor. Up to now, different dark energy models have been 
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investigated from the viewpoint of statefinder diagnos- 
tic. These models have different evolutionary trajecto- 
ries in {s, r} plane, therefore the statefinder tool can 
discriminate these models. The well known A — CDM 
model corresponds to the fixed po i nt \s = 0,r = 1} 
in the s — r plane (jSahni et al I |2003| ). The dis- 
tance of the current value of statefinder pair {s,r} 
of a given dark energy model from the fixed point 
{s ~ 0,r = l}isa valuable criterion to examine 
of model. Also the recant investigation by observa- 
tional data resulted the best fit value for statefinder 
in flat universe as {sphs = ~0-006, T-phs = 1.02} 
|M. Malekiani fc Khodam-Mohammadi 1 120123 ) . The 



corrected Friedmann equation for zero-zero component 
of field equation in flat geometry can be obtained as 
follows 



Pd) 



H 



G 
G 



(4) 



Since the value of G/G is small particularly in the late 
time accelerated universe, therefore we have ignored the 
higher time derivative of G (i.e., G/G) and also larger 
powers than one (i.e., {G/G)"^, ...). 
The equation ^ for standard model with time vary- 
ing gravitational constant can also be obtained from 
Branse- Dicke gravity by a ssuming (w = and 



other dark energy models which have been studied from j/^^^)) equation (2) of (iBaneriee fc Pavon I I2007D . 

the viewpoint of statefinder diagnostic are : ^ jjerc w is the Branse-Dicke parameter and V is Branse- 

the qu intessence DE model (ISahni et al 112003c I Alam et a l 1 ri;^i,(, gcalar field 

2003b|) , the interact ing quintessence models (|Z imdahl fc Paff )^Jconsider the derivative of G according to In a the 
2004HZhang ||2005aD. the holographic dark energy mod- ^^^ove G-corrected Friedman equation can be re-written 



els IZhang Il2005bt IZhang et al Il2007t) , th e holographic 



dark energy model in non-fl at universe (ISetare et al 



20071 ) ■ t he phantom model (|Chang et alM2007n . the 
tachyon (IShao fc Gui 1120071). the gen eralized chaply- 
gin gas model ( Malekjani et al 2011al ). the interacting 
new agegraphic DE model in flat and non-flat universe 
( Zhang iboiol : iKhodam-Mohammadi fc Malekiani Iboiol ) ,| 
the agegraphic dark energy model with and without 
interaction in flat and non-flat universe ( Wei fc Gai I 



G' 



SirGit) 



(Pm+Pd), 



(5) 



where prime is derivative with respect to a; = In a. 



2007t iMalekiani fc Khodam-Mohammadi I boiol). the 



Assuming the event horizon as an IR cut-off as 

dt f H 

-da. 



new h olographic dark energy model jMalckiani et al 



(6) 



as 



Pd 



3c2 



2011bjh . the interacting polytropic gas model (M alekiani fc KSfed&M-iMbl^gaaAadjy DE model in Eq.O is written 
2012a[l and the interacting ghost dark energy model 
(|M. Malekiani fc Khodam-Mohammadi Il2012bi ). 
The paper is organized as follows: In section (2) the G- 
corrected HDE model has been presented in fait FRW 
universe and the equation of sate Wd, deceleration pa- 
rameter q and statefinder pair {s, r} have been calcu- 
lated in the presence of time variation of G. In section 
(3) we present the numerical results and in section (4) 
the paper is concluded. 



87rGit)Rl 

In terms of dimensionless energy densities 



Pc 



Q _Pd 

^Id — , 

Pc 



(7) 



(8) 



where thepc = sj^G(t) critical energy density, the 

corrected Friedman equation([5]) can be written as 



2 The G-corrected HDE model in a FRW 
cosmology 

The Hilbcrt-Einstein action with time varying gravita- 
tional constant, G{t) — Go0(t), is 



1 



G 
G 



(9) 



5* = 



167rGo 



]d^X 



(3) 



Here we assume the scalar function (j)(t) for time depen- 
dency of G{t). Also Go is usual gravitational constant 
and Lm is the lagrangian of matter field. By varia- 
tion of above action with respect to metric g^^i, the first 



this equation is look like to the Friedman equation in 
the non-flat universe : flm + ^d = 1 — f^fe. Based on ob- 
servational data described in introduction we consider 
the negative values for ^. Therefore the added term 
G/G in right hand side of ([9]) can be interpreted as 
non-flatness parameter ilk in non-flat universe. 
In addition the evolution of Hubble parameter in terms 
of scale factor in G-corrcctcd flat universe including 
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dark matter and dark energy can be calculated from 
Eq.(|4]) as follows 

H\l - = H^[^l„,a-^ + f}rfa-3(i+--)], (10) 

where is the present value of Hubble parameter. 

The conservation equations for dark matter and dark 
energy are given by: 



+ iHp^ = 



(11) 



(12) 



Taking the time derivative of ([7]) by using Rh — 1+HRh 
and ([7]) in relation ([T^ we obtain the equation of state 
for G-corrected HDE model as follows 



Wd 



1 _ 2 vn^ 

3 3c 



1 G 
3G 



(13) 



Also, taking the derivative of (fT3|) with respect to x 
In a, we obtain 



-1 



(14) 



Here we have ignored the terms including (G/G)'^ and 

(G/G)^ and also G/G. In what follows which we derive 
and calculate, we keep only the first- order correction 
of G (i.e., G/G). 

2 

Now, derivative of 51^ = ^ — j^^i yields the 
evolutionary equation for dark energy density for G- 
corrcctcd HDE model as follows 



c E 

= -2nd\ \ — T + 11 



(15) 



In addition taking the time derivative of corrected 
Friedman equation ^ obtains 



H , 1 G , 3 , „ . 
^(l-2^) = -2(l + -<^^^) 



,G 

'G 



(16) 



Therefore the equation of motion for energy density of 
G-corrected HDE. i.e., Eq. ([T5|) is written as 



(17) 



nd i^'^) = ^dx 

(3{l + Wdnd) + ^{2-^)-3^-2 



The deceleration parameter g = — 1 — H/H^ which rep- 
resents the decelerated or accelerated phase of the ex- 
pansion of the universe, by using (|13p and (|16p . is writ- 
ten for G-corrected HDE model as 



.(l-^g)4(l + 3..^.)-^g 



(18) 



For completeness, we now derive the statefinder pair 
{s, r} in G-corrected HDE model. For this aim, by 
time derivative of (1161). we first obtain 



H 

H3 



(1---) = 
^ 2G^ 



(19) 



:(i 



:(i 



G 



Wd^di^ 



3G 



G 



4G 

){Wd^d + ^dWd) 



11 G 

Tg 



Inserting ^ and (HH) inr = -^ + 3;^ + lwe obtain 
the following equation for the parameter r of statefinder 
pair 



(1 



3G 
2G 



(20) 



9.1 ^^ M 3G, 7G\ 

-(1 + Wdnd)[wdnd{l - --) - --) 

-^(1 - ^)Wd + ^dWd) + ^§ + 1 

From ([5]), by using ([TS]) and ([20]) wc also obtained the 
parameter s in G-corrected HDE model as follows 



3/, ^ ^/ 5G, 7G\ 

-(1 + wdnd)[wdnd{i - --) ---) 



(21) 



1 



(1 



3G 
2G 



){wd^d + ^dWd) 



^G 



/ 



r3 



Wd^di^ 



3G 
2G 



5G 
4G 



In the limiting case of time-independent gravitational 
constant G (i.e., G = 0) all the above relations reduce 
to those obtained f or original ho lographic dark energy 
(OHDE) model in |Zhang Iboosi) . 



3 Numerical result 

There are many astronomical observations which show 
the time dependency of Newtonian gravitational con- 
stant. All these data are in agreement with Dyson idea 
who pointed out that G varies in the length of cos- 
mic age H^^. Based on the observational data from 
WMAP five-year observations the presen t value of Hub 



ble parameter is Ho = 6.64x10 ^^yr ^ ([Bennett et al 
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2009t IZhang fc Wu lliooj) . Moreover the astronomical 
observations are in th e line of Dirac's th eory in which 
G{t) oc i.-^ as follows |Cetto et al Ill986l ) 



G{t)^k^H{t)^k2[H{t)]-^p{ty 



(22) 



where ki and k2 are constant. The data gathered from 
SNelA data yields the be st rang for variation of G as: 



< g < l|Gaztanaga et al 1 12002[ ) and 



the data obtained from Binary Pulsar PSR1913 deter- 
mines the range of variation of g as: —(1.10 ± 1.07) x 



IQ-iiyr-^ < g < iDamour et allfim 



The data 

obtained from Helio-sesmiological determines the b est 
1.6 X 10"i2„ 



range 



yr 



-1 < g < dGuenther Ill998l) 



Another estimation for g has been done through 



astro-seismological data obtained from pulsating white 
dwarf star which yields the best range of variation as 



-2.5xlO-iV 



-1 < g < -h4.5xl0- 



yr (jBenvenuto 



20041 ). In iBiesiada fc Malec I l2004h . the range of g 



was determined as g < +4.1 x 10 ^^yr ^. It should 
be noted that all the above range of g are calculated 
for z < 3.5. Finally from the observational data of 
Big Bang nuclei-synthe sis, we have —4.0 x 10~^^yr~^ < 



g < +3 X 10 



-13„ 



yr 1 (jCopi et al 1120041) . More details 
for the variation of G with cosmic time can be seen in 



( Ray fc Mukhopadhvav 1120071 ). In previous section we 
calculated the effect of variation of G on the HDE model 
in terms of g . Therefore, we change the time derivative 



to derivative with respect to x 



In a as g 



Tjg where 



g is a dimensionless numerical value, because the di- 
mensions of Hubble Parameter is same as g . Here we 



G 
G- 



In this work wc 



1 < g < 

G ■ We 

based 



call this numerical value as a 
use the SNcIa observational data— lO^^^y?- 
which covers the other observational range of ^ 
also use the present value Hq = 6.64 x 10~^^yr 
on ob servational data from WMAP five-year o bserva- 
tions iBennett et al 1 120091: IZhang fc Wul l2009l ) . The 
parameter a, using by these observational data can be 
obtained as \a\ ^ 0.10. Therefore we choose the illus- 
trative values a = —0.1,0,0.1 which are in the order 
of the observational value. At follows we calculate the 
evolution of cosmological quantities: EoS parameter, 
energy density, deceleration parameter and statfinder 
pair of G-corrected HDE model and obtain the effect 
of parameter a on the evolution of these cosmological 
quantities. 



in Fig.(l). Here we fix the holographic parameter 
c = 0.87. Note that for this value the original HDE 
model without G correction can not enter the phantom 
regime. The black solid curve relates to original HDE 
model without G correction. The red- dashed curve 
is indicated for a = 0.1 and blue- dotted- dashed line 
represents a = 0.1. Here we see that the G-corrected 
HDE model can enter to phantom regime when a < 0, 
i.e. blue-dashed line. Hence one can conclude that the 
G-corrected HDE model can cross the phantom divide 
without a need of interaction between dark matter and 
dark energy. Also, the G-corrected HDE model crosses 
that phantom line (wd = —1) from up (wd > — 1) 
to below {wd < —1). This behavior of G-corrected 
HDE model is in agreement with recent observations 
in which the universe transits from quintessence regime 
{wd > — 1 ) to the phantorn regime {wd < —1) at the 
near past (jAlam et al 1120041 ). 



-0.7 




Fig. 1 The evolution of EoS parameter of G-corrected 
HDE model versus redshift parameter z for different illus- 
trative values of a as indicated in legend. 



3.1 EoS parameter 



By solving ([T3|) . we show the evolution of EoS pa- 
rameter of G-corrected HDE as a function of redshift 
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3.2 energy density 

Here we calculate the evolution of energy density of 
G-corrccted HDE model as a function of redshift pa- 
rameter from the early time up to late time by solv- 
ing equation ((TS|) . In Fig. (2), we plot the evolution 
of energy density fid versus of redshift for different il- 
lustrative values of a. We see that at the early times 
rid and at the late times fi^ 1, meaning the 
dark energy dominated universe at the late time. In 
this figure by fixing c ~ 0.87 the parameter a is var- 
ied as illustrative values —0.1, 0.0, +0.1 corresponding 
to dotted-dashed -blue, solid -black and dashed -red 
curves , respectively. We see that in the past times the 
dark energy becomes more dominant for positive values 
of a and at the late times the dark energy dominated 
universe can be achieved sooner for negative values. 



a' 0.5 



\ 

\\ 
■\> 

v 
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c=0.87 




a=+0.1 
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a=0 




a=-0.1 
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Fig. 2 The evolution of density parameter of Dark energy 
of G-corrected HDE model(fid) versus redshift parameter z 
for different illustrative values of a. We can see the differ- 
ent value of a result the different evolutionary trajectory in 
terms of redshift 



3.3 deceleration parameter 

Here we study the expansion phase of the universe by 
calculating the evolution of deceleration parameter q in 
G-corrected HDE model. By solving equation ((T5)) and 
using (|15p . we plot the evolution of q versus redshift 
parameter z in Fig. (3). We see that the parameter q 
start from q ~ 0.50, representing the CDM model at 
the early time. Then the parameter q becomes negative, 
representing the accelerated expansion phase of the uni- 
verse at recent epochs. Therefore the G-corrected HDE 
model can interpret the decelerated phase of the expan- 
sion of the universe at the early times and accelerated 
phase later, we fix the parameter c = 0.8 and for the 
different illustrative value of the a = —0.1, 0.0, +0.1 
corresponding to dotted-dashed -blue, sohd -black and 
dashed -red curves , respectively. We see for negative 
value of a, the accelerated expansion can be achieved 
sooner than the original HDE model(a = 0.0) and also 
positive value of a. 




Fig. 3 The evolution of deceleration parameter (g) of G- 
corrected HDE model as a function of redshift parameter z 
for different illustrative values of a. We can see accelerated 
expansion can be achieved sooner for a = +0.1 
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3.4 statefinder diagnosis 

The statefinder pair {s, r} for G-corrected HDE model 
is given by relations ([20)) and (HJ). In statefinder plane, 
the horizontal axis is defined by the parameter s and 
vertical axis by the parameter r. In Fig. (4), by putting 
nil), dm) and (Ulllin ^ and (HH and solving them 
,we obtain the evolutionary trajectories of G-correeted 
HDE model in s — r plane for different values of parame- 
ter a. By expanding the universe, the evolutionary tra- 
jectories evolve from right to left. The parameter r de- 
creases then increases, while the parameter s decreases 
forever. The trajectories cross the A—CDM fixed point 
{s ~ 0,r = 1} at the near past. In the other words, 
the G-corrected HDE model has mimicking the AGDM 
model at the near past. The present values of the 
cosmographic parameters, introduced in introduction, 
have been observationally const rained using the Markov 
Chain Monte Carlo method in ( Capozziello et aflboil ) 
as follows: Hq = 0.718, go = -0.64, tq = 1.02 , 
fco =?0.39, I = 4.05 . Using qo = -0.64 and vq = 1.02, 
we calculate the present value of statefinder parame- 
ter s as So = —0.006. Hence the observational point 
So = — 0.006, ro = 1.02 in s-r diagram is very close to 
ACDM fixed point so = 0,ro = 1. The observational 
point is indicated by green star in Fig. (4). Here we 
fix the holographic parameter c = 0.87 and vary a as 
a = —0.1, 0.0, +0.1 corresponding to dotted-dashed - 
blue, solid -black and dashed -red curves, respectively. 
We sec that different values of a result different trajec- 
tories in s — r plane. Therefore the G-corrected HDE 
model in s — r plane is discriminated for different val- 
ues of a. The colored circles on the curves represent 
the today's value of statefinder parameters {so, ^o} of 
the model. We also see that for positive values of a, 
the distance of {so,ro} from the observational point 
{sobs = —0.006, robs = 1.02} is shorter and for negative 
values of a and longer for positive values of a compare 
with original HDE model. 



4 conclusion 

In summary, we extended the holographic dark en- 
ergy (HDE) model by assuming the time depen- 
dency of Newtonian gravitational constant, G, in stan- 
dard model of cosmology. Here we obtained the G- 
corrccted Friedman equation in flat universe. Re- 
garding, the astronomical dat a from type la Su- 



pernova (IGaztanaga et al 1 120021 ) .Hulse-Tavlor Binary 
( Damour et al I Il988l). astro-seismological data from 



c=0.85 




0.7 
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Fig. 4 The he evolutionary trajectories of G-corrected 
HDE model in s — r plane for different values of a. we can see 
the different value of a result different evolutionary trajecto- 
ries, also we can see for the a = -1-0.1 the distance of present 
value from A — CDM fixed point and {sota, roi,s}(present 
value) are shorter 



2004; 



pulsa ting white dwarf stars (IBenvenuto 
20041) ■ helio-seismological data (jGuenther 



Biesiada & 
1998t) and 



the Big Bang Nuclei-synthesis data ijCopi et all 120041) . 
we obtained the parameter \a\ = ^ ~ 0.10. The evo- 
lution of EoS parameter, deceleration parameter and 
energy density parameter of HDE model in the pres- 
ence of G correction have been calculated. We showed 
that the G correction can affect the evolution of above 
parameters at the present time up to near past and 
is negligible at the early times. It was shown that 
for an illustrative value of holographic parameter c in 
which the original HDE model can not cross the phan- 
tom line, the G- corrected HDE model can achieve the 
phantom regime and cross the phantom line from up 
{wd > — 1) to below (w^ — 1) in agreement with recent 
observations (jAlam et al 1120041) . The parameter a can 
also influence on the transition from decelerated expan- 
sion to the accelerated expansion. We showed that for 
a > the transition from q > to q < earlier and for 
a < later compare with original HDE model. Finally 
we performed the statefinder diagnosis analysis for G- 
corrected HDE model and showed that the G correction 
can affect on the evolutionary trajectories of the model 
in s — r plane. We concluded that for a > 0, the dis- 
tance of present value {so,ro} from the observational 
point is shorter and for a < is longer compare with 
original HDE model, 
j^^^ T^knowledgment 
We are grateful to A. Khodam-Mohammadi for helpful 
discussions. 
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